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1. Introduction 

Muon capture on a proton is a valuable source of information about gp, the pseudoscalar 
coupling constant of the nucleon weak current One can study two processes, ordinary 
muon capture (OMC) and radiative muon capture (RMC):Q 

fi-+p-^iy^ + n (OMC), (1) 

fi'+p^u^ + n + -f (RMC). (2) 

Heavy-baryon chiral perturbation theory (HBChPT) is well suited for describing these 
processes, which involve small enough energy-momentum transfers to render the HBChPT 
series rapidly convergent. We report here our HBChPT calculation for OMC up to next- 
to-next-to-leading order. To explain the motivation and significance of our work, we first 
describe briefly the current status of OMC and RMC. 

The OMC rate in hquid hydrogen, Fuq, has been measured with 5% accuracy 0: 

r;;7 = 460±20 [sec-^]. (3) 

Theoretically, one first calculates the atomic OMC rates for a proton, Fq and Fi, where 
the suffix "0" ("1") refers to the singlet (triplet) hyperfine state of the hydrogen atom. To 
compare with experiment, one needs a theoretical framework to relate Fq and Fi to Fuq. 
For convenience, we refer to this framework as the "atom-liquid" translation formulae. A 
great deal of experimental and theoretical effort has been invested on these translation 
formulae 0. In Refs. [|l], Fq and Fi were calculated using the phenomenologically 
parametrized weak nucleon form factors with the PCAC value for gp. If we combine these 
estimates with the existing atom-liquid translation formulae, the resulting Ffj^^"^ agrees 
with F[jq^ within the experimental error 0]. In HBChPT, Fq and Fi were calculated up 
to next-to- leading order (NLO) by Bernard et al. 0]; it was also reported that a one- loop 
level calculation reproduces the known analytical PCAC correction for (/p, see Refs. P, 0. 
On the other hand, the precise empirical determination of gp is hampered by the 5% error 
in F|jq^; the sensitivity of Fj^q^ to gp is rather modest because the fixed momentum transfer 
in OMC (g^ = — 0.88m^) suppresses the contribution of the gp term, which contains the 
pion-pole structure ~ -r^— j. 

Since RMC is free from this kinematic constraint, it can be a more sensitive probe of gp, 
despite its extremely small branching ratio. A recent TRIUMF experiment succeeded in 
measuring dVYmc/ dE^, the absolute photon spectrum for RMC in liquid hydrogen [§, H. 
If one uses atomic RMC amplitudes calculated in the phenomenological minimal coupling 



method ||T0[, and if one adopts the existing atom-liquid translation formulae, then the 
observed dVjiyiQ/ dE^ cannot be reproduced unless gp is artificially increased from its 
PCAC value by as much as 50 %. However, such a large deviation of gp from its PCAC 
value is extremely unlikely according to an HBChPT calculation H . This very astonishing 



feature reported for dV-^yiQ/ dE^ motivated reexamination of the formalism used in |]T0 



to calculate the RMC amplitude, and several calculations based on HBChPT have been 



In this article OMC and RMC always refer to capture in a hydrogen target. 
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carried out [|, [TT], [T^. A next-to-leading order (NLO) calculation [|TI| indicates that 



HBChPT essentially reproduces dVnuc/dE^ given in |T^. A next-to- next-to-leading order 
(NNLO) calculation ||T^ has confirmed that the HBChPT expansion converges rapidly and 
that loop corrections to dT^Mc/dE^ are tiny. Furthermore, a recent calculation |^ that 
incorporates the explicit A degrees of freedom into a tree-diagram HBChPT calculation 
suggests that the inclusion of the A modifies the spectrum only by 5%, a result consistent 
with the earlier finding of Beder and Fearing [0. This change is not large enough to 



remove the above-mentioned anomalous gp value. Thus the systematic analyses based 
on HBChPT strongly indicate that no drastic changes in the atomic RMC amplitudes 
from the existing estimates should be expected. It then seems likely that the dTjiMc/dE^ 
problem is caused by the currently adopted atom-liquid translation formulae. 

Meanwhile, an experiment that uses a hydrogen gas target to directly measure Pq with 
1% accuracy is planned at PSI [|T^. The use of the gas target eliminates ambiguities due 
to the molecular capture processes. The envisaged 1% accuracy will significantly increase 
precision with which the empirical value of gp is determined. We note, however, that, to 
make comparison between theory and experiment at the 1% level, the existing estimate 
of Pq based on an HBChPT calculation up to NLO needs to be improved. First, one 
must ascertain that the input physical constants like /^r and gA have sufficient precision. 
Secondly, NNLO loop corrections need to be evaluated^. 

In this article we present an HBChPT calculation for OMC up to NNLO in which the 
influence of uncertainties in the low-energy constants is carefully examined. In addition to 
Po and Pi, we calculate the neutron helicity asymmetry (to be defined later). It is found 
that Pq obtained here is significantly larger than the previous estimates [0, ^, whereas 
Pi essentially agrees with the literature values. The larger value of Pq is primarily due 
to the updated larger value of gA- A second main point of our paper is to discuss the 
observational ramifications of our new estimates. If we use the new values of Pq and 
Pi together with the "standard" atom-liquid translation formulae [H], the resulting rf^^°^ 



turns out to be significantly larger than T'^^J' in Eq. (1^), another possible indication 



liq 

liq 

that the existing atom-liquid translation formulae may require reexamination. We shall 
argue that this difficulty is probably related to the above-mentioned "anomaly" in RMC, 
and that there is possibility to resolve these two problems simply and simultaneously by 
invoking a molecular mixing parameter discussed by Weinberg 111 



2. Heavy-baryon chiral perturbation theory (HBChPT) 

HBChPT is a low energy effective field theory of QCD, which has a systematic per- 
turbative expansion in powers of Q/A^, where Q is a small typical four-momentum 
scale characterizing a process in question, or the pion mass m^r; A^^, is the chiral scale, 
A^ ~ 47r/^ ~ rriAT ~ 1 GeV. A typical scale Q in muon capture (both OMC and RMC) is 
the muon mass = 105.7 MeV, and hence Q/A^ ~ 0.1. One therefore expects a rapid 
convergence of relevant chiral perturbation series for muon capture; the explicit HBChPT 
calculations 0, 0, [1^, |1^, |18[ are consistent with this expectation. 



^ The transition amplitude of OMC was calculated in HBChPT up to NNLO by Fearing et al. [|7|, 
but the capture rate was not explicitly given in that work. 
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The effective chiral Lagrangian is expanded as 

C = Y.C, = C^ + C^ + C2 + ---. (4) 

The subscript v denotes the order of terms, z/ = (i + n/2 — 2, where n is the number of 
nucleon hues and d the number of derivatives or powers of involved in a vertex. The 
terms relevant to our calculation are 

£0 = N[iv-D + 2igAS-/S\N + fli:i[-/\-/\ + ^y (5) 

A = ^N[{v-Df-D^ + 2gA{v-^,S-D}-i{l + h,)[S^,S^]f^p]N, (6) 



where is the l/mj^ Lagrangian given in Ref. ||T2|. Furthermore 



D, = d, + T^, (8) 
A. = lW,d,^}-'-eF^^+'-^Fi:e, (10) 



= ex^^±^X^^, X = ml e = exp(zf-7r/2A). (12) 

In these expressions ■u^ and are the isovector vector and axial-vector external fields, 
respectively; = (1, 0) is the velocity four- vector, and = (0, cr/2) is the nucleon spin 
operator. We ignore the isospin breaking effect and use rriN = (mp + mn)/2 as the nucleon 
massQ. Our effective Lagrangian contains the low energy constants (LECs), 65, C3, C13 and 
Ci4.[| The LECs, 65, C13 and C14, are finite constants fixed by experiments. To one-loop 
order, 65 = Ky + g\m.„m]\i/4:7if^, where =3.706 is the isovector anomalous magnetic 
moment. The constant C13 is fixed by the mean square axial radius deduced from (anti) 
neutrino-proton scattering, (r^i) = 0.65 fm |l^; its numerical value is C13 = (47r/^)^-^^-|^ 
= 4.88. The parameter C14 is fixed by the GT discrepancy defined by 

. _ 9AmN _ 2ml ^ 



^ In calculating the transition amplitudes we ignore isospin breaking but, in evaluating the phase 
space, we do retain the neutron-proton mass difference (see later). 

^ Our notations for the LECs are different from those in Ref. Ml . The relations between them are 



&5 = C6, C3 = Sio, Ci3 = 2(477/^)^^24/5^, Ci4 = B; 



23- 
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where g^^N is the ttNN couphng constant. The one-loop diagrams are renormahzed by C3. 
Integrating in d = 4 — 2e dimension, we have 



C3 



1 + 5g\)R + cf , with R 



1 ^ , 47r/i2 

- + 1 - 7 + In- 



(14) 



where 7 = 0.5772 ■ ■ -, and the mass scale /i is a parameter in dimensional regularization. 
Note that we include /i into c§ to avoid the n dependence in the amplitudes. The param- 
eter is fixed by the empirical radius of the isovector Dirac form factor (r}^)^ = 0.585 



fm^ lEOi. Thus, from 



£2_ 



;i.267). 



967r2/2 



we deduce 



5.39 (5.42) with c/a = 1-26 



3. Atomic OMC rates and neutron helicity asymmetry: Formalism 

The OMC process is effectively described by the current-current interaction. Thus the 
transition amplitude reads 



M 



(15) 



where G^jVud = Gp is the Fermi constant, is the leptonic weak current, and J" is 
the nucleon weak current. The leptonic current is simply given by = 'Uj/7a(l — 75)^/^, 
whereas J" is a much more complex object refiecting hadron dynamics. Here we evaluate 
J° in HBChPT up to NNLO (one-loop) chiral order. 

Since in HBChPT the nucleon current J° is expanded in terms of I/itln, it is convenient 
to write J" in the Pauli-spinor form. The time and spatial components of the nucleon 
current J" = Jy — are written as 



4iQ) = fiil)^ Mq) = ^'^ X qf^iq) + qf^iq), 

A{q) = ^ ■ qffiq), Uq) = ^fiiq) + q{^ ■ q)f2{q), 



(16) 
(17) 



where we have suppressed the initial- and final-nucleon spinors as well as the common 
factor 2miy. The form factors, and {i = 1,2,3), introduced here may be referred 
to as the non-relativistic (NR) form factors. The relations between the NR form factors 
and the standard relativistic nucleon weak form factors are given in the Appendix. Note 
that f'Y and are in fact related by current conservation. 

The NR form factors calculated up to one-loop order in HBChPT read 



V 



1 + 

4 



4 o 1 + I7gl 



;q 



q' + 



1 



(AnUy^ 18(471^)2-^ Am% 
' ^lil + 2gl)ml-'-±^q^ 



-i^ + ^v]q 



V 



(47r/.)2 
1 



/o(g), 



2m N 

k1 



1 + Ky) + 



q' + 



9A 



n{Aml - q^) 
4m^ 



mo(g) 



2m7v ' 



(18) 
(19) 
(20) 
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9a 
9a 

9A 



1 + 



Cl3 



2(47r^) 



N , 



2m 



2(47rA)2 

- fl - A.(g)g2) |g1 



5m 



N 



(21) 

(22) 
(23) 



where A^(g) ^ = (q'^ — m^) ^ is the renormahzed pion propagator, and the one- loop 
functions are given by 



mo{q) 



dx ln[l — x(l — x)g^/m^], 



dx 



(24) 
(25) 



'1 + x(l — x)(p /ml. 

The total OMC rate from a muonic hydrogen atom in a hyperfine state S is given as 



r5 = i^-^;^El^(^;^'^^ 



2 

I 5 



with 



K = 



(26) 



(27) 



The helicity amplitude M(/i; S", 5*2) in Eq. (^) is specified by the final neutron helicity h 
[h = L, R), the initial hyperfine spin S, and its 2;-component, Sz- In Eq. (^), Ei, is the 
final neutrino energy given by^ 



E. 



(m^ + mp)^ - m^ 
2(mu + m„) 



99.15 [MeV]. 



(28) 



The factor 0^(0) appearing in K is the value at the origin of the radial wavefunction for 
the fi~-p ground state; thus |0^(O)p = (amp^/vr, where a is the fine structure constant 
and m^ = m^rrip/ (m^ + rUp). 

When the neutron helicity is monitored, we define the neutron helicity asymmetry as 



as 



Ts{L) - Ts{R) 
Vs{L) + Vs{R) 



(29) 



where Ts{L) {Vs{.R)) is the rate of OMC from an initial atomic spin state, 5", leading to 
a final left-handed (right-handed) neutron.^ 

In calculating the capture rates for different neutron helicities, it is convenient to choose 
the direction of the emitted neutrino as the z-axis. In general there are eight helicity 



^See footnote!^. 

^ Thus, if the neutron hehcity is not monitored, Ts — ^s{L) + Ts{R)- 
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/.^ 




LO 


1.000 








1.267 


-0.426 





NLO 





0.248 


0.053 








0.045 


NNLO 


-0.030 


-0.004 





-0.021 


0.006 





Total 


0.970 


0.244 


0.053 


1.246 


-0.419 


0.045 



Table 1: Numerical values of the OMC form factors in Eqs. ([16|) and ([T7|), calculated for 
each chiral order with the use of qa = 1.267 and g.j,N = 13.4. 



amplitudes, but with this particular choice we have only three non-zero amplitudes. They 
are 

M(L; 0, 0) = ^(/r + 2f^ + + 3/i^ + fi + f^), (30) 

M{L- 1, 0) = - + - fi + f2 + f^l (31) 

M{R- 1, -1) = + - ft - ft - fi), (32) 

where [3 = AGpmN\Jrn^j,Ey. Correspondingly, we have Tq{L) = |M(L; 0, 0)^, ro(-R) = 
0,ri(L) = fs:||M(L;l,0)p, and Vi{R) = Kl\M{R;l, 

4. Atomic OMC rates and neutron helicity asymmetry: numerical results 

As emphasized above, at the level of precision of our concern, we need to be particularly 
careful about the accuracy of the input physical parameters. The most updated value of 
Gi3 is Gp/V2 = 0.8030 ± 0.0008 x 10"^ GeV-^. For qa and Q-kN, we use as the standard 
values gA = 1.267 ||21| and g-^N = 13.4. 

Since the momentum transfer for OMC is fixed, we calculate the NR form factors for 
\q\ = Ey and = — 2m^Ey. The results for each order of HBChPT expansion are 
given in Table |l|. The table clearly shows that the chiral perturbation series converges 
very rapidly. As for the helicity-dependent amplitudes, we obtain 

M(L; 0, 0) = 3.45/5, M(L; 1, 0) = -0.77/3, M(i?; 1, -1) = 0.15 /3. (33) 



In Table |] we give our numerical results for the atomic capture rates, Fq and Fi, 



along with HBCHPT calculations of Bernard et al. The second column labeled 

"This work" represents the results obtained in our NNLO-calculation. Comparing this 
column with the third column that gives the results of our NLO-calculation, we note that 
the NNLO corrections decrease Fq significantly (3.9%). Thus it is clear that, in order to 
achieve theoretical precision that matches the 1% accuracy of the planned PSI experiment 

HI 



on Fo 



one must take into account the NNLO terms. 



We now turn to comparison with the other calculations quoted in Table The esti- 
mates of PrimakofT ||l| and Opat are based on the phenomenological parameterization 
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This work 


This work 


Bernard et a/.|^ 


1 Bernard et al.^ 


Primakoff [0] 


Opat § 




(NNLO) 


(NLO) 


(NNLO) 


(NLO) 




To 


695 


722 


687.4 


711 


664±20 


634 


Ti 


11.9 


12.2 


12.9 


14.0 


11.9±0.7 


13.3 



Table 2: Comparison of calculated atomic OMC rates. Tq (Fi) is the capture rate [in 
sec~^] for the initial singlet (triplet) hyperfine state. The entries for the columns labeled 
"This work (NNLO)" and "This work (NLO)" have been obtained with the use of qa = 
1.267 and g^N = 13.4. 



of the weak nucleon current with gp fixed at its PCAC value; Primakoff retained the 
relativistic kinematics, whereas Opat used a non-relativistic expansion of the amplitudes 
in terms of l/mjy- Bernard et a/.'s results come from an NLO-calculation |^ and an 
NNLO-calculation ||2^. The results of the two NNLO-calculations, the second and fourth 
columns in Table exhibit a 1 % difference. Bernard et al. used a different value of gA 
{gA = 1.26). They also expanded the phase space and atomic wavefunction in Eq. ( p7D 
in powers of m^/mjv. For the sake of comparison, suppose that, in our NNLO calculation 
we apply expansion of m^/mjsi to the phase space, while keeping the exact atomic wave- 
function, but ignoring neutron-proton {n-p) mass difference. Then Fq calculated up to 

l/m^ 

1/m^ would become Fg ^ = 708 sec~^, which is to be compared with Fq = 695 sec~^ in 
Table ^. If we expand both phase space and wavefunction, ignoring n-p mass difference, 
we obtain r^'"'^ = 705 sec^^. The different l/m^v expansion schemes show only a tiny 
numerical difference, but, the rates in both cases are enhanced by ~ 2 % compared with 
the exact phase space case. This is due to ignoring the nucleon mass difference; the rate 
is proportional to and Ey increases by 1 % when the n-p mass difference is ignored. In 
the present work we retain the exact final phase space and atomic wave-function;P^ viz., 
we apply HBChPT expansion only to the transition amplitudes. On the other hand, the 
1 % difference between the two NNLO-calculations in Table |^ may stem from differences 
in input parameters as well. 

According to Table ^ Fq's calculated in HBChPT (our present result and those of 
Bernard et al.) have significantly larger values than the earlier theoretical estimates. This 
is primarily due to the fact that the modern HBChPT calculations employ an updated 
value of gAi which is larger than the older values. Primakoff |jl[ used gA = 1.24, while 
Opat [Q used = 1-22. In order to illustrate the sensitivity of our results to the input 
physical parameters, we show in Table ^the values of Fq and Fi corresponding to different 
values of gA and g-^N] the LEC Ci4, which is determined by gA and g-^N via Eq. ([T3|) , is 
also listed. As can be seen from the last four rows in the table, for a given value of 
gAi variations in (jf^jy causes only minor changes in the capture rates; even though these 

^'^ We do not, however, include the finite-nucleon-size effect on the atomic wave- function; this effect is 
known to reduce the capture rate by ~0.4 %, see e.g., Appendix f in Ref. |^ . 
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9a 


9ttN 


Ci4 


To 


Ti 




1.22 


13.40 


-2.59 


656 


11.3 


0.865 


1.24 


13.40 


-2.07 


672 


11.6 


0.893 


1.26 


13.40 


-1.54 


689 


11.9 


0.918 


1.26 


13.05 


-0.65 


692 


11.6 


0.927 


1.267 


13.40 


-1.36 


695 


11.9 


0.925 


1.267 


13.05 


-0.47 


698 


11.7 


0.934 



Table 3: The LEG Ci4, atomic capture rates Tq and Fi [sec and the neutron hehcity 
asymmetry ai, calculated in HBChPT up to NNLO for various choices of qa and g^^N. 



variations lead to a difference of a factor of ~3 in Ci4, the corresponding changes in the 
rates in the last four rows are modest; 1.3% for Tq and 2.5% for Fi. Thus the most crucial 
input parameter here is qa- 

Finally we discuss the helicity asymmetry. Due to the V — A weak interaction the final 
neutron is purely left-handed when an initial atomic state is in the hyperfine singlet state 
(5 = 0); thus we have ao = 1 as a trivial identity. For the initial hyperfine triplet state 
{S = 1), the final neutron can have both left- and right-handed helicity components, and 
therefore ai can have a non-trivial value; our calculation gives ai = 0.925. Thus almost 
all out-going neutrons are polarized left-handedly, a result consistent with Weinberg's 



observation 1 5 



5. OMC rates for liquid hydrogen targets 

As mentioned, in order to relate theoretical estimates of the atomic OMC rates to the 
capture rate measured in liquid hydrogen, one needs (in our terminology) the "atom- 
liquid" translation formulae. We briefly describe here the atomic and molecular physics 
input that underlies these formulae. A muon stopped in liquid hydrogen quickly forms a 
muonic atom (/i-p) in its Bohr orbit. The atomic hyperfine triplet state (S=l) is rapidly 
transformed into the singlet state (S=0); this hyperfine transition rate is known to be 
Aio ~ 1.7 X 10^° sec~^ indeed a very large value. A muonic atom and a hydrogen 
molecule collide with each other to form a p-^-p molecule, predominantly in its ortho 
state (with the two proton spins parallel to each other). Let App^ be the rate of transition 
from the atomic hyperfine singlet state to the ortho p-fi-p molecular state. Meanwhile, 
the ortho p-fj.-p molecular state decays to the lower-lying para p-fi-p molecular state. Let 
Xop stand for this decay rate. Taking into account these atomic and molecular processes, 
one relates Fuq, the OMC capture rate in liquid hydrogen, to the atomic capture rates 
(Fo and Fi) via the formula 



liq 



^" ro+ . . ^\ (rom + T,^^], (34) 



Ao + App^ Aq + App^ Ao + Aop \ /\o / 

where Aq is the muon decay rate, Aq = 0.455 x 10^ sec~^. In this equation Tom (Xpm) 
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represents the rate of muon capture from the ortho (para) p-fi-p molecular state. These 
rates are usually calculated using the formula, 

^on^ = (^Fo + ill) , Tpm = 2^p + ^Fi) . (35) 

The factors 270 and 27P account for modifications of the muon wavefunction as it changes 
from the atomic Bohr orbit to the p-^-p molecular state; according to Ref. 0, 270 = 1.009 
and 27p = 1.143. The validity of Eq. (p^) will be discussed later in the text. 



As for App^, there are several conflicting experimental results, see e.g. Ref.[Tj]. We 
quote here the lowest and highest reported values: A^^^ = 1.89 ± 0.20 x 10^ sec~^, and 
\IZ = 2.75 ± 0.25 X 10^ sec-^ The current theoretical estimate A*'!,, = 1-8 x 10^ sec"^ 
lies near the lower edge of the lowest experimental value, which is a rather uncomfortable 
situation. Furthermore, the current experimental and theoretical values for Aop do not 
agree with each other: = 4.1 ± 1.4 x 10^ sec'^ El and A*^ = 7.1 ± 1.2 x 10^ sec'^ 



0. Thus it seems fair to say that the existing "atom-liquid" formulae are not totally free 
from uncertainties and that these ambiguities can affect our interpretation of the OMC 
rate in liquid hydrogen. 

In the following we use for App^ the value adopted in Ref. App^ = 2.5 x 10^ sec^^; 
for Aop, we use two representative values: A*p and A^p^. Using the OMC rates of Primakoff 
in Table |] and A*^, Bakalov et al. obtained Fuq = 490 ± 10 sec~^, in good agreement 
with the data, V^^^ = 460 ± 20 sec~^. However, if we use the values of Fq and Fi obtained 

th 



m our NNLO-HBChPT calculation together with A*^, then Eq. (||) gives a much larger 



op' 

rate, Fuq = 518 sec~^. If we adopt Aop = Aop^, Tug becomes even larger: Tug = 532 sec^^. 
Thus, the use of the updated values of Fq and Fi as obtained here in combination with 
the commonly used "atom-liquid translation formulae" spoils the previously reported good 
agreement between T\l^^"^ and Fjj^^. So, in addition to the problem of the RMC photon 
spectrum discussed in the Introduction, another serious problem seems to be lurking in 
the OMC sector. 

6. A mixing molecular parameter to fit the OMC and RMC data 

In view of the fact that these two problems occur in the experiments involving liquid 
hydrogen targets, it seems of interest and of importance to reexamine the reliability of the 
formulae hitherto used in the literature to relate the atomic capture rates to Fuq. Although 
a thorough investigation of this issue is beyond the scope of this article, we wish to discuss 
here a particular aspect of molecular physics input which seems relevant to the present 
issue but so far has not been fully examined. Taking up an early observation made by 



Weinberg [TH] , consider the possibility that in liquid hydrogen two ortho molecular p-fx-p 
spin states, S'p^p = 1/2 and 3/2 may be populated. If this indeed happens, Fo^ in Eq. 
( P^D should be replaced with 

rL = erom(i/2) + (1 - Oro,n(3/2), (36) 

where Fom(l/2) = Fom of Eq. ( ^5l) and Fom(3/2) = 270F1. According to Weinberg [jlSl , 



the mixing parameter ^ can be in the range of 0.5 < ^ < 1. Only theoretical estimates of ^ 
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1.0 0.95 0.90 0.85 0.80 


riiq(Aop) 
riiq(A^^^) 


518 499 480 461 442 
532 512 492 472 452 



Table 4: The OMC rate in liquid hydrogen, Fuq [sec calculated for five different values 



of the molecular mixing parameter ^, Eq. (^6]), and for two typical choices (explained in 
the text) of Xop, the ortho-to-para transition rate in a p-^-p molecule. The second (third) 
row corresponds to the use of A*p = 7.1 xlO"^ sec 



iK7 



4.1 xlO^ sec" 



exist, and its literature value is ^ ^ 1 [0, |[; to our knowledge ^ has never been measured 



experimentally. To study the consequences of ,^ 7^ 1, we have calculated Fuq with the use 
of Eqs. (|3^ ) ,(|35| ) ,(|36|), for several values of ^. For Fq and Fi we have used the results 
of our NNLO calculation while for Xop we have considered the two representative values 
discussed above: A*p and A^^^. The results are shown in Table ^. The table indicates that, 
if 0.8 < ^ < 0.9, the theoretical and experimental values of F^q are in good agreement. 
We next argue that the introduction of ^ in this range leads to resolution of the RMC 
problem as well. We first remark that Eqs. (^4D,(|35D,(|36D, can be used, mutatis mutan- 
dis, for RMC as well, in particular, for calculation of dT^Mc/dE^, the photon spectrum 
for RMC in liquid hydrogen. With the atomic RMC transition amplitudes previously 
obtained in a NNLO-HBChPT calculation |T2[, we have evaluated dTnuc/dE^ for vari- 



ous values of ^; the other atomic and molecular population parameters are kept fixed at 
the values used in Ref. 0]: 6.1% atomic hyperfine singlet state, 85.4% ortho p-fi-p state 
and 8.5% para p-fi-p state. The results are shown in Fig. ||. The dashed line (lowest 
curve) represents the no-mixing case, ^ = 1, which corresponds to dT^Mc/dE^ obtained 
in |]12|. For comparison , we also show in the figure (the solid line) the result obtained 



in a modified version of the Fearing model [TQ] wherein gp is taken to be 1.5 times the 



PCAC value. This line represents the best fit curve to the observed dTjiMc/dE^ in the 
analysis reported in p|, One can see from the figure that ^ in the range of 0.8-0.9 leads 
to a photon spectrum that is satisfactorily close to the "observed" spectrum (solid line) 
for E^ > 60 MeV. 

7. Discussion and Conclusions 

We have considered in our HBChPT calculation up to NNLO contributions. Some 
remarks on possible higher order effects are in order here. As discussed in Ref. a one- 
loop diagram in N'^LO which contains a vertex with the anomalous magnetic moment, 
KviTLfj,/mN ~ 0.5, can be comparable in size to the NNLO diagrams. This means that the 
correction to the capture rate due to the N^LO terms may reach the 1-2 % level. Thus, for 
a more precise theoretical prediction of the OMC rates, it can in principle be important 



If one replaces Eq. ( p4[ ) with an older and simpler "atom-liquid" formula of Ref. |23|, the OMC 
data can be fit with ^ = 1, and the RMC data with ^ ~ 0.95 within the framework of an NLO-HBChPT 



calculation, see Ref. 26 
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Figure 1: The RMC photon spectrum, dTRuc/dE^. The sohd hue labeled "1.5PCAC" 
is the result one would obtain from Fearing's model if the value of gp is assumed to be 
1.5 times the standard PCAC value. The other curves represent the results of NNLO- 
HBChPT calculations with ^ = 1.0, 0.9 and 0.8. 



to include the N^LO corrections, although the existing uncertainty in the values of qa 
and Qt^n may not warrant the effort. Moreover, at the level of N^LO, the isospin breaking 
effects as well as QED corrections are expected to give sizable contributions. 



In conclusion, we have carried out a HBChPT calculation of the atomic OMC rates 
to next-to-next-to- leading order. Our result indicates that, once the measurement of the 
hyperfine-singlet atomic OMC rate reaches 1% accuracy, as envisaged in a PSI experi- 
ment, theoretical predictions based on HBChPT must include at least NNLO corrections. 
Furthermore, we have shown that both the OMC rate and the RMC photon energy spec- 
trum measured in liquid hydrogen can be reproduced by introducing the molecular mixing 
parameter ^ in the range of ^ ~ 0.8 to 0.9. It seems interesting to examine whether this 
range of ^ is realistic. 

Finally, we have shown that the neutron helicity asymmetry for OMC from a hyperfine 
triplet state is ~93%. 
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APPENDIX 

Assuming the absence of the second-class current, one can express the nucleon vector and 
axial currents in terms of four form factors: 



J, 



Un{p'] 



G'y(g)7" + GM(g) 



2m 



■N 



m, 



(37) 
(38) 



where Gy(g), Gm{(i), GA^q), and Gp{q) are the vector, weak-magnetism, axial-vector, 
and pseudo-scalar form factor, respectively; q = p' — p and Up (m„) is the Dirac spinor for 
the proton (neutron). These standard relativistic form factors are related to the NR form 
factors, /^'s and ff's, defined in Eq. Up to 0{l/m%) 



= Gv{q) 1 - 



Gniq), 



2m 



N 



[Gviq) + GMiq)]\q\, = 



1 



N . 



2m 



■N 



Gp{q) I 



Gv{q)\q 

.2 



2mnm]\[ 



1 + 



^m% 



'GA(q) + ^^^^]\q 
V ^ 2m^mjsi ) 
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